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Abstract 

We consider the elliptic-elliptic, focussing Davey-Stewartson equations, which have an explicit bright line soliton 
solution. The existence of a family of periodic solitons, which have the profile of the line soliton in the longitudinal 
spatial direction and are periodic in the transverse spatial direction, is established using dynamical systems 
arguments. We also show that the line soliton is linearly unstable with respect to perturbations in the transverse 
direction. 

Resume 

Solitons periodiques du systeme de Davey-Stewartson elliptique-elliptique focalisant. Nous considerons 
les equations de Davey-Stewartson focalisantes dans le cas elliptique-elliptique, lorsqu’elles possedent une solu¬ 
tion unidimensionnelle de type soliton. En utilisant des methodes de la theorie des systemes dynamiques, nous 
montrons I’existence d’une famille de solutions bidimensionnelles qui out le profil d’un soliton dans la direction 
spatiale longitudinale et sont periodiques dans la direction spatiale transverse. Nous montrons egalement que le 
soliton unidimensionnel est lineairement instable vis-a-vis des perturbations transverses. 


1. Introduction 

The Davey-Stewartson equations 

iAt + SAxX + Ayy + {1 i\A\^ + D2(t>x)A — 0, (1) 

l3(j>xx + (t>yy - 13\A\1 = 0, ( 2 ) 

where £ = ±1, 71, 72, 73 G IR\{0} with 72+73 = ±2, arise in the modelling of wave packets on the surface 
of a three-dimensional body of water; the variables A = A{x, y, t) and (j) = ^(a:, y, t) are the complex wave 
amplitude and real mean flow and the signs of the parameters depend upon the particular physical 
regime under consideration (see Ablowitz & Segur [2, §2.2]). In the literature the cases 71 + 72 = 2 
and 71+72 = —2 are termed respectively focussing and defocussing, and the system is classified as 
hyperbolic-hyperbolic, hyperbolic-elliptic, elliptic-hyperbolic or elliptic-elliptic according to the signs of e 
and 73. Certain special cases of the mixed-type systems are often referred to as DS-I and DS-II and are 
known to be completely integrable (see Ablowitz & Clarkson [1, p. 60]). Note that solutions of (1), (2) 
which are spatially homogeneous in the j/-direction satisfy the cubic nonlinear Schrodinger equation 

iAt + Axx + (71 + 72)|^P^ = 0 (3) 

(where f is recovered from (2)). 
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Solutions of (3) which converge to an equilibrium as a; —)■ ±oo and are 27r-periodic in t are referred 
to as line solitons. In the defocussing case the equation admits a ‘dark’ line soliton which decays to a 
nontrivial equilibrium, while in the focussing case it has the ‘bright’ line soliton 

t) = e’* sech(x) (4) 

which satisfies A*{x,t) —>■ 0 as a; —)■ ±oo. In this note we examine periodic solitons which decay as 
X — )■ ±oo and are periodic in y and t, and in particular consider how they emerge from line solitons 
in a dimension-breaking bifurcation. Explicit formulae for dark periodic solitons have been obtained for 
the integrable versions of the equations by Watanabe & Tajiri [9] and Arai, Takeuchi & Tajiri [3]; here 
we establish the existence of bright periodic solitons to the elliptic-elliptic, focussing equations (e = 1, 
7 i + 72 = 2, 73 > 0) under the additional condition 72 > 0 by dynamical-systems methods. 

Theorem 1.1 Suppose that £ = 1 , 71 -I- 72 = 2 and 72, 73 > 0 . There exist an open neighbourhood Af of 
the origin in M, a positive real number loq and a family of periodic solitons {e^*Us{x,y),(j)six,y)}seM' lo 
(1), (2) which emerges from the bright line soliton in a dimension-breaking bifurcation. Here 

Us{x,y) =sech(a;) -i-u'^{x,y), (j)s{x,y) = tanh(a;) + (j)'^{x,y), 

in which u^(- , •), (/)^(- , •) are real, have amplitude 0(|s|) and are even in both arguments and periodic in 
their second with frequency ujq -I- 0(|sp). 

We also present a corollary to this result which asserts that the bright line soliton is transversely linearly 
unstable and thus confirms the prediction made by Ablowitz & Segur [2, §3.2], 

Theorem 1.2 Suppose that £ = 1 , 71 -I-72 = 2 and 72, 73 > 0 . For each sufficiently small positive value 
of X the linearisation of (1), (2) at A*{x,t) = e'‘sech(a;), (j)*{x,y) = tanh(a;) has a solution of the form 
e^*'^'^{A{x,y),(f>{x,y)), where {A{x,y),(j){x,y)) is periodic in y and satisfies {A{x,y),(j){x,y)) —>■ ( 0 , 0 ) as 
X —>■ ± 00 . 

In the remainder of this article we suppose that £ = I, 71 -I- 72 = 2 and 72, 73 > 0. Equations (I), (2) 
with these coefficients arise when modelling water waves with weak surface tension. The existence and 
transverse linear instability of periodic solitons for the water-wave problem in this physical regime has 
recently been established by Groves, Sun & Wahlen [6]. 

2. Spatial dynamics 

The equations for solutions of (1), (2) for which A{x,y,f) = e'*(^ui{x,y,t) -\-iu 2 {x,y,t)) (and u\, U 2 are 
real-valued) can be formulated as the evolutionary system 


Uly=Vl, (5) 

Vly = U 2 t - Uixx +Ui- (71 -I- 71M2 + l 2 (l>x)ui, (6) 

U2y = V2, (7) 

V2y = -Uit - U2xx +U2- (71 Ml 4- 71^2 + 724>x)u2, ( 8 ) 

fy = f), (9) 

i’y = -ISffxx + + u\)x, ( 10 ) 


where the spatial direction y plays the role of time. To identify an appropriate functional-analytic setting 
for these equations, let us first specialise to stationary solutions, so that 
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Uly=VX, ( 11 ) 

Viy = -Uixx + Ml - (71M1 + 71U2 + l 2 (px)ui, (12) 

U2y = M2, (13) 

V2y = -U2xx +U2- ( 71^1 + 71^2 + l24>x)u2, (14) 

<Py = '>P, (15) 

i’y = -734>xx+73iuj+ul)x. (16) 

Equations (11)-(16) constitute a semilinear evolutionary system in the phase space X = x 


L^(IR) X X L^(IR) x x L^(IR); the domain of the linear part of the vector field defined 

by their right-hand side is H = iJ^(IR) x x x iJ^(IR) x _ff^(IR) x This evolution¬ 
ary system is reversible, that is invariant under y i—)■ —y, (mi, mi, M2, M2, V') > 5 '(mi, vi, M2, W2, <(', V’)) 

where the reverser S : X ^ X is defined by S'(mi, mi, M 2 , M 2 , V') = (mi,—mi, M 2 ,—M 2 , </>,— ■*/')• K is 

also invariant under the reflection R : X ^ X given by R{ui{x),Vi{x),U2{x),V2{x), (l>{x),'ip{x)) = 
(ui(—x),vi(—x),U2(—x),V2(—x), —(/>(—x), —^^(—x)), and one may seek solutions which are invariant un¬ 
der this symmetry by replacing X and £> by respectively 

Xr :=XnFixR = iJi(M) x L^(R) x x L^(R) x ff^(R) x L^(R) 

and 

D, := BnFixR = II^(R) x ff^(R) x II^(R) x ff^(R) x ff^(R) x ff^(R), 

where 


ffg(R) = {w € H^(R) : w(x) = w(—x) for all x G K}, 
ff^(R) = {w€ Ef”(M) : w(x) = -w(-x) for all x G R}. 

It is also possible to replace Dy by the extended function space 

:= H!{R) X iJi(M) X H!{R) x H^(R) x J/,%(M) x H^{R), 

where 

H'^^^{R) = {im G ifo(.(M) : Wx G H^{R.),w{x) = —w{—x) for all x G M} 

(a Banach space with norm ||m ;||*^2 := ||M>a;||i)- This feature allows one to consider solutions to (11)-(16) 
whose (((-component is not evanescent; in particular solutions corresponding to line solitons fall into this 
category (see below). 

Each point in phase space corresponds to a function on the real line which decays as x —>■ 00 , and 
the dynamics of equations (11)-(16) in y describes the behaviour of their solutions in the y-direction. In 
particular, equilibria correspond to line solitons (the equilibrium 

(mi(x), Mi(x), M 2 (x), M 2 (x), (j)*{x), ■0*(x)) = (sech(x), 0,0,0, tanh(x), 0) 

corresponds to the line soliton (4)), while periodic orbits correspond to periodic solitons (see Figure 1). 
In Section 4 we construct dimension-breaking bifurcations by writing 

(Mi,Mi,M2,M2,<(),f/') = {Ul,Vl,<l>*,U2,V2,(j)*,tp*) + (m) , m) , M 2 , M^ ,(()', V’O (17) 

and seeking small-amplitude periodic solutions of the resulting evolutionary system 

Wy = Lw + N{w) (18) 

for w = (m), m(, M 2 , M 2 , (/>', V’O) where 
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Figure 1. A family of periodic solutions surrounding a nontrivial equilibrium solution to (11)-(16) in its phase space (left) 
corresponds to a dimension-breaking bifurcation of a branch of periodic solitons from a line soliton (right, plot of u{x,y)). 
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and we have dropped the primes for notational simplicity. Note that (18) has the invariant subspace 
X = {(m 2 ) M 2 ) = (0, 0)}, and we define = X^Ci X, n X and l)* = U* (7 X. 

Returning to (5)-(10), observe that these equations constitute a reversible evolutionary equation with 
phase space H^{{—to,to),X); the domain of its vector field is H‘^{{—to,to),X) 7 H^{{—to,to),D) and its 
reverser is given by the pointwise extension of S' : X X to H^{{—to,to),X). Seeking solutions of the 
form (17), we find that 

Wy = Twt + Lw + N{w), 


where T(ui, ui, U 2 , M 2 , (/i,'0) = (0, M 2 ,0, —mi, 0, 0) and we have again dropped the primes. In Section 5 
we demonstrate that the solution (u{, u*, uj, mJ, 0*, 0*) of (5)-(10) is transversely linearly unstable by 
constructing a solution of the linear equation 


Wy = Twt-\-Lw (19) 

of the form e^*u\{y), where u\ S Cj((M,X) 7 Ch{^,D) is periodic, for each sufficiently small positive 
value of A. 










3. Spectral theory 


In this section we determine the purely imaginary spectrum of the linear operator L : D C X ^ X. To 
this end we study the resolvent equations 

{L — ikl)w = (20) 

for L, where w = (ui, vi, U 2 , U 2 , <(','*/’), and A: € M \ {0}; since L is real and 

anticommutes with the reverser S it suffices to examine non-negative values of k, real values of ui, U 2 , 
cj), vl, and purely imaginary values of u\, u\, ui, U 2 , V'- Observe that (20) is equivalent to the 
decoupled equations 

v\ + iku\ 

-|- ifc(()^ 

X L2(M) ^ L2 


where Ai : 
defined by 



{A2 + k'^I)u2 = v\+ iku^, 

X L2(M) and A 2 ■■ i?2(ffi) C L2( 




are 


Ai 


-Uixx +Ui- (371 -I- 72) sech^(a;)ui - 72 sech(a:)^^ 
-l 3 <Pxx + 273 (sech(a;)ui) 2 : 

^2^2 = —U2xx +U 2 — 2sech^(x)M2; 


the values of vi, V2 and ip are recovered from the formulae 

vi = u\+ ikul, V2 = u\+ iku2, Ip = (p'^ + \k(p. 

It follows that L — iki is (semi-)Fredholm if Ai + k"^! and A 2 + k^I are (semi-)Fredholm and the dimension 
of the (generalised) kernel of L — \kl is the sum of those of Ai + k^ I and A 2 + k'^I. 

Lemmata 3.1 and 3.2 below record the spectra of Ai and A 2 ; part (i) of the following proposition (see 
Drazin [4, Chapter 4.11]) is used in the proof of the former while the latter follows directly from part (ii). 

Proposition 3.1 

(i) The spectrum of the self-adjoint operator 1 — df — 6sech^(a;) : H'^fR.) C L'^fR) —>• L'^(R) consists of 
essential spectrum [l,oo) and two simple eigenvalues at —3 and 0 (with corresponding eigenvectors 
sech^(a;) and secii (x)). 

(ii) The spectrum of the self-adjoint operator 1 — df — 2sech^(a;) : H'^(R) C L‘^(R) —>• L'^(R) consists of 
essential spectrum [l,oo) and a simple eigenvalue at 0 (with corresponding eigenvector sech(a;)J. 
Lemma 3.1 The spectrum of the operator Ai consists of essential spectrum [0,oo) and an algebraically 
simple negative eigenvalue — Wg whose eigenspace lies in Ll{R) x 

Proof. First note that Ai is a compact perturbation of the constant-coefficient operator 
H^{R) X H^{R) C L^{R) X L^{R) L^{R) x L^{R) defined by 

iui,(p) (-Ml 

XX + Ui, -J 3 (p XX ) ; 

whose essential spectrum is clearly [0,oo); it follows that (Jess(.4i) = [0, 00 ) (see Kato, [8, Chapter 
IV, Theorem 5.26]). Because Ai is self-adjoint with respect to the inner product {{ul,(p^),{u‘l,(p‘^)) = 
{ul,ul)o + 5727 ^^(^^, for L^(IR) x L^(IR) the remainder of its spectrum consists of negative real 
eigenvalues with finite multiplicity. 

One finds by an explicit calculation that 

2 T2 f 2 

{Ai{ui,(p),{ui,<p)) = {ui-uixx-Qsech {x)ui,ui)o-\-— / {(px - 2 sech{x)ui) dx, 

^ Jr 
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which quantity is positive for {ui,(f>) G x with 

{{ui,(j)), (sech^(a;),0)) = 0 

(see Proposition 3.1(ii)). It follows that any subspace of x upon which Ai is strictly negative 

definite is one-dimensional. The calculation 

lim (^i(sech(x), 2(()i^(a;) tanh(x)), (sech(x), 2(()i^(a;) tanh(a:))) = — 
fl—>00 3 

where (j){R) = x(x/i?) and y € is a smooth cut-off function eqnal to unity in [—1,1], shows that 

inf (t{Ai) < 0, so that the spectral subspace of x H^{R) corresponding to the part of the spectrum 

of Ai in (—oo, —e) is nontrivial and hence one-dimensional for every sufficiently small value of e > 0. We 
conclude that Ai has precisely one simple negative eigenvalue —loq. 

Finally, the same argument shows that |l 2 (r)xL 2 (r) also has precisely one simple negative eigen¬ 
value. It follows that this eigenvalue is — Wg, whose eigenspace therefore lies in Te(K) x Lq(IR). □ 
Lemma 3.2 The spectrum of the operator A 2 consists of essential spectrum [1, 00) and an algebraically 
simple negative eigenvalue at 0 whose eigenspace lies in Lf. (M). 

Corollary 3.3 The purely imaginary number ik belongs to the resolvent set of L for fc G K \ {0, iwg} 
and iiwg are algebraically simple purely imaginary eigenvalues of L whose eigenspace lies in X^. 


4. Application of the Lyapunov-Iooss theorem 

Our existence theory for periodic solitons is based upon an application of the following version of the 
Lyapunov centre theorem for reversible systems (see looss [7]) which allows for a violation of the classical 
nonresonance condition at the origin due to the presence of essential spectrum there (a feature typical of 
spatial dynamics formulations for problems in unbounded domains) provided that the ‘looss condition at 
the origin’ (hypothesis (viii)) is satisfied. 

Theorem 4.1 (looss-Lyapunov) Consider the differential equation 

Wr = L{w) + N{w), (21) 

in which w{t) belongs to a real Banach space X. Suppose that y, Z are further real Banach spaces with 
the properties that 

(i) Z is continuously embedded in y and continuously and densely embedded in X, 

(a) L : Z C X ^ X is a closed linear operator, 

(in) there is an open neighbourhood U of the origin in y such that L G L{y,X) and N G Cg^(7/, A) 
(and hence N G ^(7/ n Z,X)) with N{0) = 0, d7V[0] = 0. 

Suppose further that 

(iv) equation (21) is reversible: there exists an involution S G T(X) n T(y) H £(Z) with SLw = —LSw 
and SN(w) = —N{Sw) for all w GU, 

and that the following spectral hypotheses are satisfied. 

(v) iiwg are nonzero simple eigenvalues of L; 

(vi) inujo G p{L) for n G Z\{ —1,0,1}; 

(vii) \\{L - muJoI)~''-\\x-^x = o{l) and \\{L - muJoI)~^\\x-^z = 0{1) as n± 00 ; 

(viii) For each GU the equation 

Lw = —N{w^) 

has a unique solution w Gy and the mapping ^ w belongs to Cg^(7/,3f). 
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Under these hypotheses there exist an open neighbourhood I of the origin in M and a 
continuously differentiable branch {(v(s),a;(s))}sg/ of reversible, 2 t:/ uj{s)-periodic solutions in 
C'pgj.(M, 3^ 0 A") n C'per(K, iV © 2) to (21) with amplitude 0(|s|). Here the direct sum refers to the decom¬ 
position of a function into its mode 0 and higher-order Fourier components, the subscript ‘per’ indicates 
a 2 t:/ uj{s)-periodic function and a;(s) = wq © 0(|sp). 

Theorem 1.1 is proved by applying the looss-Lyapunov theorem to (18), taking X = X„, y = £)*, Z = 
D„ and U = (and of course t = y and S(ui,vi,4>,il)) = {ui, —vi, —'tp)). The spectral hypotheses 
(v) and (vi) follow from Corollary 3.3, while (vi) and (vii) are verified in respectively Lemma 4.2 and 4.3 
below. 

Lemma 4.2 The operator L satisfies the resolvent estimates \\{L — ikl)~^\\^ = 0(|A:|“^) and 

\\{L - = 0(1) as \k\ oo. 

Proof. Notice that 

L{ui,vi,(l),'ip) = (Bi{ui,vi),B2{(l>,'ip)) -\-C{ui,vi,(l),f^), 
where Bi{ui,vi) = (wi,+ Ui), ^ 2 (<(',V’) = and 

C(mi,Ui,(/),V') = ( 0,-(371 0 72 )sech^(a;)ui - 72 sech(a:)(/)a;, 0,-73()) 0 273(sech(x)'Ui)^). 

Writing = Xi x X 2 , = Yi x Y 2 and equipping Xi with the usual inner product and X 2 with 

the inner product ((^^, ^^), ((/)^, ^/>^)) = observe that Bj : Yj C Xj —>• Xj is 

self-adjoint, so that 

\\{B,-ikI)-^\\^^^^^<\k\-^ 
for fc 0. Furthermore \\Bj{-)\\^, = || • Wy., so that 

\\iB,-ikI)-^^^^y^ = \\B,iB,-ikI)-^^^^j,^ = \\I + ikI{B,-ikI)-^^^^j,^ < 2 

for fc 0. It follows that B = Bi x B 2 ■ Dr ff —>■ X„ satisfies the estimates 

\\{B-ikI)-^y^^j,^<\k\-\ \\{B-ikI)-^y^^f,^<2 ( 22 ) 

for fc 0 . 

Finally, note that C : Xr —> is bounded, whence 

\\CiB-ikI)-^^^^^^=0{\k\-^) 

as |A:| —>■ 00 . Consequently I — C{B — ikI)~^ ■. Xy. ^ Xy. is invertible for sufficiently large values of |A:| with 

11(7 - C{B - ifc/)-i)-i= 0(1) (23) 

as |A:| —>■ 00 , and the stated result follows from the identity 

{L-ikl)-^ = {B-ikI)-^{l-C{B-ikI)-^)~^ 

and the estimates (22), (23). □ 

Lemma 4.3 The equation 

Lw = -N{w^) (24) 

has a unique solution w G O* for each ■u;! € O* and the formula i-G w defines a smooth mapping 

O* —> 0*. 

Proof. Equation (24) is equivalent to the equations 

ui - uixx - 6 sech^(cc)Mi = /(w'^'), 
where /(w'^) = (371 0 72 ) sech(a;)(M|)^ 0 ^ 2 u\(j>\ 0 71 (u|)^, and 

= (wj)^ 0 2 sech(x)Mi, Ui = 0 , Ip = 0. 

The result thus follows from Proposition 3.1(i) and the fact that / and (■ui,m|) ^ (wj)^ 0 2sech(a;)ui 
dehne smooth mappings O* —>■ Ll{R) and Hf{R) x H)(R) —)• □ 
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5. Transverse linear instability 


Finally, we demonstrate the transverse linear instability of the line soliton using the following general 
result due to Godey [5]. 

Theorem 5.1 (Godey) Consider the differential equation 

VT = Tvt + Lv, (25) 

in which v(t, t) belongs to a real Banach space X. Suppose that y, Z are further real Banach spaces with 
the properties that 

(i) L : Z C X ^ X and T -.ycX^X are closed linear operators with Z cy, 

(a) the equation is reversible: there exists an involution S G Z{X) n L{y) n L{Z) with LSv = —SLv 
and TSv = —STv for all v G Z, 

(Hi) L has a pair Ziloq of isolated purely imaginary eigenvalues with odd algebraic multiplicity. 

Under these hypotheses equation (25) has a solution of the form e^^vx^r), wherev\ G (^^(M, T) nC(M, Z) 
is periodic, for each sufficiently small positive value of X; its period tends to 2 t:Iu}q os A —> 0. 

Theorem 1.2 is proved by applying Godey’s theorem to (19), taking X = X,y = X and Z = D (and of 
course t = y and S(ui,vi,U2,V2,4>,4’) = (ui,—vi,U2,—V2,4>,—'f)). The spectral hypothesis (hi) follows 
from Gorollary 3.3. 
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